DFTT 18/2007 
KIAS-P07045 

Zooming in on light relic neutralinos by direct detection and measurements of galactic 

antimatter 

A. Bottino,i'B F. Donato,i'Q N. Fornengo,!'! and S. Scopel^'l 

^ Dipartimento di Fisica Teorica, Universita di Torino 
Istituto Nazionale di Fisica Nucleare, Sezione di Torino 
via P. Giuria 1, 1-10125 Torino, Italy 
^Korea Institute for Advanced Study 
Seoul 130-722, Korea 
(Dated: February 2, 2008) 

The DAMA Collaboration has recently analyzed its data of the extensive WIMP direct search 
(DAMA/Nal) which detected an annual modulation, by taking into account the channelling effect 
which occurs when an ion traverses a detector with a crystalline structure. Among possible implica- 
tions, this Collaboration has considered the case of a coherent WIMP-nucleus interaction and then 
derived the form of the annual modulation region in the plane of the WIMP-nucleon cross section 
versus the WIMP mass, using a specific modelling for the channelling effect. In the present paper 
we first show that light neutralinos fit the annual modulation region also when channelling is taken 
into account. To discuss the connection with indirect signals consisting in galactic antimatter, in 
our analysis we pick up a specific galactic model, the cored isothermal-sphere. In this scheme we 
determine the sets of supersymmetric models selected by the annual modulation regions and then 
prove that these sets are compatible with the available data on galactic antiprotons. We comment 
on implications when other galactic distribution functions are employed. Finally, we show that 
future measurements on galactic antiprotons and antideuterons will be able to shed further light on 
the populations of light neutralinos singled out by the annual modulation data. 

PACS numbers: 95.35.-|-d,11.30.Pb,12.60. Jv,95.30.Cq 



I. INTRODUCTION 



In a recent paper [l[ the DAMA Collaboration has 
analyzed the data of its extensive WIMP direct search 
(DAMA/Nal) ^ which measured an annual modulation 
effect at 6.3 a C. L., by taking into account the chan- 
nelling effect. This effect occurs when an ion traverses 
a detector with a crystalline structure In Ref. [l| 
implications of channelling have been discussed in terms 
of a specific modelling of this effect for the case of the 
DAMA Nal(Tl) detector; it is shown that the occur- 
rence of channelling makes the response of this detec- 
tor to WIMP-nucleus interactions more sensitive than in 
the case in which channelling is not included. Therefore, 
when applied to a WIMP with a coherent interaction 
with nuclei, the inclusion of the channelling effect im- 
plies that the annual modulation region, in the plane of 
the WIMP-nucleon cross section versus the WIMP mass, 
is considerably modified as compared to the one derived 
without including channelling. The extent of the modi- 
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fication depends on the specific model-dependent proce- 
dure employed in the evaluation of the channeling effect 

Q. 

These properties are shown in Fig. [1] where the quan- 
tity CT^"aiai°'^ denotes the WIMP-nucleon scalar cross- 
section, ^ ~ pwiMp/po is the WIMP local fractional 
matter density and is the WIMP mass. The dashed 
line denotes the annual modulation region derived by the 
DAMA Collaboration without including the channeling 
effect 0. The solid line shows the annual modulation 
region derived by the same Collaboration when the chan- 
neling effect is included as explained in Ref. [H . 

The regions displayed in Fig. [T]arc derived by varying 
the WIMP galactic distribution function (DF) over the 
set considered in Ref. 0| and by takinjg into account other 
uncertainties of different origins [11, [J. Fig. □ shows that 
the effect of taking channelling into account is that the 
annual modulation region modifies its contour with an 
extension towards lighter WIMP masses. Most remark- 
ably, for WIMP masses < 30 GcV, the WIMP-nucleon 
cross section involved in the annual modulation effect 
decreases sizeably, up to more than an order of mag- 
nitude. As mentioned before, the specific shape of the 
annual modulation region depends on the way in which 
channelling is modelled 

These features are of great importance for a specific 
dark matter candidate, the light ncutralino, which was 
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FIG. 1: WIMP-nucleon scattering cross-section as a function 
of the WIMP mass. The solid (dashed) hne denotes the an- 
nual modulation region derived by the DAMA Collaboration 
with (without) the inclusion of the channeling effect. The two 
regions contain points where the likelihood- function values 
differ more than 4a" from the null hypothesis (absence of mod- 
ulation). These regions are obtained by varying the WIMP 
galactic distribution function (DF) over the set considered 
in Ref. [J and by taking into account other uncertainties of 
different origins [l| . The scatter plot represents supersymmet- 
ric configurations calculated with the supersymmetric model 
summarized in the Appendix. The (red) crosses denote config- 
urations with a neutralino relic abundance which matches the 
WMAP cold dark matter amount (0.092 < ^x^^ < 0.124), 
while the (blue) dots refer to configurations where the neu- 
tralino is subdominant {Q,^h^ < 0.092). 



extensively investigated in Refs. 0,0, Si- Actually, in 
these papers we analyzed light neutralinos, i.e. ncutrali- 
nos with a mass < 50 GcV, which arise naturally in 
supersymmetric models where gaugino mass parameters 
are not related by a GUT-scale unification condition. In 
Refs. 0,01 it is proved that, when R-parity conserva- 
tion is assumed, these neutralinos are of great relevance 
for the DAMA/Nal annual modulation effect. In these 
papers it is also shown that in MSSM without gaugino 
mass unification the lower limit of the neutralino mass is 
m^>7 GeV 

In Fig. [ll superimposed to the annual modulation re- 
gions is the scatter plot of the supersymmetric configu- 
rations of our model, whose features are summarized in 
the Appendix. One sees that, also when the channeling 
effect is taken into account, the light neutralinos of our 
supersymmetric model fit quite well the annual modula- 
tion region. 

In the present paper we consider the phcnonicnological 



consequences for light neutralinos when the annual mod- 
ulation region is the one indicated by the solid line in 
Figd] More specifically we examine the properties of our 
supersymmetric population of light relic neutralinos in 
terms of the possible antimatter components generated 
by their pair annihilation in the galactic halo. 

To do this, we have to resort to a specific form for the 
WIMP DF. We take as our representative DF a standard 
cored isothermal-sphere model, though we do not mean 
to associate to this model prominent physical motivations 
as compared to other forms of DFs. Analyses similar to 
the one we present here for the cored isothermal-sphere 
can be developed for other galactic models. We will com- 
ment about some of them, selected among those consid- 
ered in Ref. (we will follow the denominations of this 
Reference to classify our DFs). 

The scheme of the present paper is the following. In 
Sect. II, we show how the model presented in Refs. 
@,0jS| fits the DAMA/Nal annual modulation regions of 
Ref. [l| for the case of the cored isothermal-sphere model. 
In Sect. Ill we combine these results with constraints 
derivable from available data on cosmic antiprotons; we 
also discuss the sensitivity of upcoming measurements on 
cosmic antiprotons for investigating the neutralino popu- 
lations selected by the annual modulation regions. Com- 
plementary investigations by measurements of galactic 
antideuterons are presented in Sect. IV. Conclusions are 
drawn in Sect.V. The main features of the supersym- 
metric scheme adopted here are summarized in the Ap- 
pendix. 



II. THE ANNUAL MODULATION REGION IN 
VARIOUS HALO MODELS 

As mentioned above, in the present paper we take 
the cored isothermal sphere as the representative model 
for our detailed evaluations. Similar analyses can be 
developed for other galactic models; we will comment 
about some of them. The density profile of the cored- 
isothermal sphere (denoted as Evans logarithmic model, 
or Al model, in Ref. 0]) is: 



p{r) 



AttG {Rl + r2)2 



(1) 



where G is the Newton's constant, vq is the local value 
of the rotational velocity and Rc is the core radius. 

The value Rc = 5 kpc will be used for the core ra- 
dius. For the parameter vq we will consider the values 
vq = 170,220,270 km sec~^, which represent the min- 
imal, central and maximal values of vq in its physical 
range 10[. For each value of vq, we will consider the 



3 




FIG. 2: WIMP-nucleon scattering cross-section as a function of the WIMP mass. The solid contours denote the DAMA/Nal 
annual modulation regions for a cored isothermal halo, derived by including the channeling effect with the model explained 
in Ref. The different panels refer to different galactic halo-model parameters, according to the analysis of Ref. 0]: «o 
is the local rotational velocity, po is the local dark matter density. The scatter plot shows the configurations for neutralino- 
nucleon scattering in gaugino non-universal supersymmetric models. The (red) crosses denote configurations with a neutralino 
relic abundance which matches the WMAP cold dark matter amount (0.092 < Q^h^ < 0.124), while the (blue) dots refer to 
configurations where the neutralino is subdominant {fl^h^ < 0.092). 



minimal and the maximal values of the local dark matter 
density, pjf™ and p™'^^, as determined in Ref. Q. Then, 
specifically, we will discuss the following sets of values: 
i) vo = 170 km sec"! with p^'"' = 0.20 GeV cm'^ and 
^max ^ 0.42 GeV cm-3; ii) Vo = 220 km sec'^ with 
p;,'^'" = 0.34 GeV cm-3 and p^'''' = 0.71 GeV cm^^; iii) 
Vo = 270 km sec'^ with p^'"^ = 0.62 GeV cm'^ and 
pgiax ^ 07 GeV cm-3. 

Now, we turn to a comparison of the DAMA/Nal an- 



nual modulation regions of Ref. |3| with the theoretical 
predictions of our supersymmetric model. Fig. [2]displays 
the DAMA annual modulation regions in the case of the 
Al model |ll| : the various insets correspond to the rep- 



resentative values of the parameters Vo and po previously 
defined. The regions of Fig. [2] are derived by the DAMA 
Collaboration from their data of Ref. taking into 
account the channelling effect and under the hypothe- 
sis that the WIMP-nucleus interaction is coherent. They 
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represent regions where the likehhood- function values dif- 
fer more than 4a from the null hypothesis (absence of 
modulation) [l^]- The scatter plot is the same as in Fig. 
[TJ It is remarkable that light ncutralinos are able to pro- 
vide a good fit to the experimental data. This occurs for 
values of vq and po which arc in the low-medium side of 
their own physical ranges, i.e. vq ~ (170 - 220) km sec^^ 
and po ~ (0.2 - 0.4) GeV cm~^. The hght neutralinos 
involved in this fit stay in the mass range ~ (7 — 30) 
GcV. We remark that in case the channelling is not in- 
cluded, the DAMA regions would be sizeably displaced 
as compared to the ones displayed in Fig. 2, similarly 
to what is shown in Fig. 1. An example of the effect 
of including channelling in the determination of the an- 
nual modulation region, when the Al model for the DF 
is taken, is explicitly displayed in Fig. 7 of Ref. [1]. As 
already remarked before, in connection with Fig. 1, for 
WIMP masses < 30 GeV, the WIMP-nucleon cross sec- 
tion involved in the annual modulation effect decreases 
sizeably, when channeling is included, up to more than 
an order of magnitude. This implies that, not including 
channelling, the fit of the experimental data with light 
neutralinos would require values of po (po — (0.6 — 1) 
GcV cm^'^) higher than the ones previously derived. Also 
Vq would be in the high side of its physical range. These 
properties are of relevance for the implications which will 
follow. 

When the channelling effect is taken into account and 



no rotation of the halo is considered, it turns out 
that the features of the annual modulation region in the 
771^ — ^0'""^^'}'™" plane do not differ much when the galactic 



DF is varied, for many of the galactic DFs considered in 
Ref. Thus, for instance, for a matter density with 
a Navarro- Frenk- White profile (A5 model of Ref. [^) 
or for an isothermal model with a non-isotropic velocity 
dispersion (Bl model of Ref. Q) the physical situation 
is very similar to the one depicted in Fig. [21 However, in 
the case of DFs with triaxial spatial distributions (within 
the class D of Ref. 0]) and for models with a corotating 
halo there can be an elongation of the annual modulation 
region towards heavier masses Further insight into 
the properties of light neutralinos are expected from the 
future results of the DAM A/LIBRA experiment [1]. 

We wish also to stress that the distribution of WIMPs 
in the Galaxy could deviate from the models mentioned 
above, mainly because of the presence of streams. For 
modification of the annual modulation region in these in- 
stances see Ref. Q • It is worth mentioning that in the nu- 
merical derivation given above, also uncertainties of other 
origin may intervene. Suffice it to mention that the size- 
able uncertainties which affect strength of the coupling 
of the neutralino to the nucleon . 

In this paper, among the searches for WIMP direct de- 



tection we discuss only the DAMA/Nal experiment, since 
this is the only experiment having at present the capa- 
bility to measure the annual modulation effect, which is 
a distinctive feature for discriminating the signal against 
the background in a WIMP direct search [l5|. For up- 
dated reviews about other experiments of WIMP direct 
detection, see Ref. 1161. 



III. GALACTIC ANTIPROTONS 

As shown in Ref. among the various searches for 
indirect signals due to annihilation of light WIMPs, the 
cosmic rays antiprotons provide the most significant con- 
straints. For this reason we now examine how this kind 
of limits applies to the light neutralinos singled out by 
the DAMA/Nal annual modulation regions. 

So-called secondary antiprotons are produced in the 
Galaxy via interaction of proton and helium cosmic rays 
with the interstellar hydrogen and helium nuclei. A thor- 
ough calculation of the secondary antiproton spectrum 
has been performed in Ref.[l7[, where the antiprotons 
generated by spallation processes are propagated using a 
two-zone diffusion model described in terms of five pa- 
rameters. Two of these parameters, Kq and i5, enter the 
expression of the diffusion coefficient: 



K = Ko(3R'^, 



(2) 



where R is the particle rigidity. The other three param- 
eters are the Alfen velocity Va, the velocity of the con- 
vective wind Vc, and the thickness L of the two large 
diffusion lay ers which sandwich the thin galactic disk. 
In Ref. 17[ it is shown that the experimental antipro- 



ton spectrum is fitted quite well by the secondary com- 
ponent from cosmic-rays spallation (with a, = 33.6 
with 32 data points), calculated with the set of the diffu- 
sion parameters which is derived from the analysis of the 
boron-to-carbon ratio (B/C) component of cosmic-rays. 
The values of this set of best-fit parameters (denoted as 
median), together with their 4cr uncertainty intervals, are 
given in Table I. The theoretical uncertainty on the dif- 
fusion parameters reflects into a (10 - 20)% uncertainty 
on the calculated spectrum of secondary antiprotons. 

Primary antiproton fluxes can be generated by anni- 
hilation of neutralino pairs. We have evaluated these 
fluxes for the supcrsymmetric configurations selected by 
the annual modulation regions, i.e. the light neutralino 
populations which stay inside the annual modulation re- 
gions displayed in the insets of Fig. [S] These correspond 
to the cases: A) vq = 170 km sec~^, po 
GeV cm^"^; B) vq = 170 km sec~^, po 
GeV cm-3; C) vq 
GeV cm-3. 
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FIG. 3: Antiproton flux at p kinetic energy Tp = 0.23 GeV, as a function of the WIMP mass and for a cored isothermal halo. 
Each raw correspond to a different set of cosmic-rays propagation parameters: the upper, median and lower rows refer to 
the set which provides the maximal, median and minimal antiproton flux, according to the analysis of Ref. (3]. The light 
gray points denote conflgurations with a neutralino-nucleon scattering cross section outside the corresponding DAMA/Nal 
allowed region. The bold (colored) points refer to configurations compatible with the DAMA/Nal regions. These last points 
are further differentiated as follows: (red) crosses denote conflgurations with a neutralino relic abundance which matches the 
WMAP cold dark matter amount (0.092 < Q^h^ < 0.124), while (blue) dots refer to configurations where the neutralino is 
subdominant {Q,xh? < 0.092). The solid horizontal line shows the maximal allowable amount of antiprotons in the BESS data 
[2o| over the secondary component; the dashed and dotted lines denote estimates of the PAMELA and AMS sensitivities to 
exotic antiprotons for 3 years missions, respectively. 



The antiproton fluxes originated in the dark halo by 
the neutralino pair-annihilation processes have then been 
propagated in the diffusive halo using the three sets of dif- 
fusion parameters (minimal, median and maximal) given 
in Table I. The procedure for the evaluation of these 
fluxes is the one illustrated in Refs. 



18, m 



As 



it was shown in Ref. [l8|, the uncertainty in the dif- 
fusion/propagation parameters, contrary to the case of 
the secondary antiprotons, induces a large uncertainty 
on the primary flux. 

To show quantitatively how the experimental data can 
constrain our supersymmetric configurations, in Fig. [3] 
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TABLE I: Astrophysical parameters of the two-zone diffu- 
sion model for galactic cosmic-rays propagation, compatible 
with B/C analysis [l3l and yielding the maximal, median and 
minimal primary antiproton flux. 



we display the (top-of-the-atmosphere) antiproton flux 
evaluated at a specific value of the antiproton kinetic en- 
ergy, Tp = 0.23 GcV, for the three populations A, B, and 
C defined above. The shaded (yellow) region denotes the 
amount of primary antiprotons which can be accommo- 
dated at Tp = 0.23 GcV without entering in conflict with 
the BESS experimental data [i^l and secondary antipro- 
tons evaluations 17| . The dashed horizontal line denotes 
our estimated sensitivity of the PAMELA detector [2l| to 
exotic antiprotons after a 3 years running: it corresponds 
to the admissible excess within the statistical experimen- 
tal uncertainty if the measured antiproton flux consists 
only in the background (secondary) component. The es- 
timate has been performed by using the background cal- 
culation of Ref. 



and refers to a l-cr statistical un- 
certainty. All the supersymmetric configurations in Fig. 
[3] above the dashed line can be potentially identified by 
PAMELA as a signal over the secondaries, while those 
which are below the dashed curve will not contribute 
enough to the total flux in order to be disentangled from 
the background. The dotted horizontal line represents a 
similar estimate, but referred to the AMS detector 
for a 3 years data-taking. Crosses (red) and dots (blue) 
denote ncutralino configurations selected by the annual 
modulations regions and with 0.092 < O^/i^ < 0.124 and 
il^h^ < 0.092, respectively. Faint (gray) dots represent 
configurations which are outside of the annual modula- 
tion regions. Fig. [3]shows that, for values of the diffusion 
parameters close to the minimal set, present antiprotons 
data do not set constraints. For values of the diffusion pa- 
rameters around the median set, we have that: in case A, 
most of the neutrahno configurations with fi^ft,^ < 0.092 
and a few with 0.092 < ft^h^ < 0.124 remain uncon- 
strained; in cases B and C only subsets of neutralinos 
with riy^h^ < 0.092 survive. When the diffusion param- 
eters approach the values of the maximal set, only very 
few SUSY configurations survive in case A. From Fig. [3] 
we see that the possibility of exploring our relevant neu- 
trahno configurations by future measurements of galac- 
tic antiprotons (PAMELA and AMS) is quite good. As 
was remarked in Sect. II, in case the channelling is not 
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FIG. 4: Areas of compatibility between the annual modu- 
lation regions of Fig. [2] and the antiproton data for a neu- 
trahno of mass of 20 GeV, plotted in the parameter space 
defined by the height of the diffusive halo L and the rigidity- 
dependence parameter S of the diffusion coefficient of Eq. (0). 
The galactic halo model is a cored-isothermal sphere. The 
dots, (red) squares and (blue) circles refer to: vq = 170 km 
sec"\ po = pjf'" = 0.20 GeV cm"^, vq = 220 km sec"\ 
Po = P™'" = 0.34 GeV cm"^, and vo = 170 km sec"\ 
po — Po^^'' = 0.42 GeV cm~^, respectively. Each set of points 
shows the region in the L-5 plane which fits at 99.5% C.L. 
the antiproton data of BESS [20| . 



included, the fit of the experimental data of annual mod- 
ulation with light neutralinos would imply values of po 
higher than the those characterizing the sets A, B and 
C, previously defined. This property, together with the 
fact that the antiproton flux depends on the square of poj 
might cause tension between the annual modulation data 
and the constraints implied by present measurements of 
galactic antiprotons. 

In the previous analysis, we have considered the an- 
tiproton flux evaluated at a specific value of the antipro- 
ton kinetic energy, Tp = 0.23 GeV. This analysis can be 
extended by examining the properties of a global fit to 
the full low-energy antiproton spectrum, using the same 
procedure which was used in Ref. 



19[. We mentioned 



above that this spectrum is fitted quite well by the sec- 
ondary component from cosmic-rays spallation 17|. As 
a conservative criterion for constraining our supersym- 
metric configurations, we can perform a analysis on 
the antiproton data for the supersymmetric configura- 
tions compatible with the annual modulation study. As 
an illustrative application of this analysis, for the super- 
symmetric configurations compatible with each set A, B 
and C, we calculate the antiproton fiux for all the prop- 
agation parameter combinations which keep the B/C fit 
within a 4cr of uncertainty [2^. The primary and sec- 
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ondary fluxes are then required to flt the antiproton data 
at 99.5% C.L. (x^ <60). The results of this calculation 
are displayed in Fig. d] for = 20 GeV in the plane 
of the two diffusion parameters L and 5. We see that, 
depending on the specific isothermal-sphcrc parameters, 
we identify different regions of compatibility of the an- 
tiproton signal in terms of the astrophysical parameters 
which govern the diffusion of galactic cosmic-rays. 



IV. GALACTIC ANTIDEUTERONS 

Formation of antideuterons in cosmic rays proceeds 
through production of an antiproton and an antineutron 
pair by spallation (secondary production) or by WIMP 
pair annihilation (primary production) [2j]. The coa- 
lescence process of antiproton and antineutron is easier 
in WIMP annihilation, since their parent particles are 
at rest in the galactic frame. Therefore, at low energies 
the primary spectrum is much enhanced as compared to 
the secondary one 24|. This feature makes the search 
for antideuterons particularly attractive for an indirect 
investigation of WIMPs ^Ell. 

The fluxes of the antideuterons produced in the dark 
halo by the neutralino pair-annihilation processes have 
been calculated following the method described in Ref. 
[23 | and propagated in the diffusive halo using the three 
sets of diffusion parameters (minimal, median and max- 
imal) given in Table I. 

In Fig. [5] we display the (top~of-the-atmosphere) an- 
tideutcron flux evaluated at the value Tp ~ 0.23 GeV/n 



for the three population A, B, and C defined above. The 
notations for the scatter plot are as in Fig. [3l that is: 
crosses (red) and dots (blue) denote neutralino configu- 
rations selected by the annual modulations regions and 
with 0.092 < Vt^K^ < 0.124 and O^/i^ < 0.092, respec- 
tively; faint (gray) dots represent configurations which 
are outside of the annual modulation regions. The hor- 
izontal dashed and dotted lines denote estimated sensi- 
tivities to antideuterons of the GAPS and AMS ^\ 
detectors. 

Fig. [5] shows that measurements of galactic an- 
tideuterons are perspectively very promising for investi- 
gating our light neutralino populations. Moreover, when 
antideuteron data will become available together with 
the antiproton ones, correlations in the two data sets 
will provide strong confidence in a possible presence of 
a signal, as can be appreciated by comparing Figs. 2 
and 3: for most of our relevant light-mass neutralinos, 
a signal should be present both in the antiproton and 
antideuteron channel. 



V. CONCLUSIONS 

In the present paper we have considered the annual 
modulation regions which the DAMA Collaboration has 
recently determined, by including also the channelling 
effect which occurs when an ion traverses a detector 
with a crystalline structure, such the detector of the 
DAMA/Nal experiment. The inclusion of the channelling 
effect implies that the annual modulation region is con- 
siderably modified as compared to the one derived with- 
out including channelling. The extent of the modification 
depends on the specific model-dependent procedure em- 
ployed in the evaluation of the channeling effect. 

In the present paper we have considered the phe- 
nomenological consequences for light neutralinos when 
the annual modulation region includes the channelling 
effect as modelled in Ref.[l|. We have proved that these 
annual modulation data are fitted by light neutralinos 
which arise naturally in supersymmetric models where 
gaugino mass parameters are not related by the a GUT- 
scale unification condition. 

The precise range of the neutralino mass which fits 
the annual modulation data depends on how the WIMP 
galactic distribution function is modelled and on a num- 
ber of other assumptions, such as those mentioned in 
Sect. II. As an example, we have worked out in detail the 
case of a cored isothermal sphere DF. For this instance, 
the neutralino mass stays in the range ~ (7 — 30) 
GeV, for values of the local rotational velocity, vq, and of 
the local dark matter density, po, in the low- medium side 
of their own physical ranges, i.e. vo — (170 - 220) km 
sec^^ and po ~ (0.2 — 0.4) GeV cm~'^. Similar ranges are 
found also in the case of a Navarro-Frcnk- White profile 
or for an isothermal model with a non-isotropic velocity 
dispersion. 

We have then shown that the populations of light neu- 
tralinos selected by the annual modulation regions are 
consistent with present data on galactic antiprotons. We 
have also derived the intervals of the diffusion param- 
eters which provide this agreement in correlation with 
the specific galactic halo model. For instance, for neu- 
tralinos with a mass of 20 GeV and a cored isothermal 
model with vq = 170 km s~^ we have 0.55 < 5 < 0.85 
and L < 3 kpc when po = pg'^'^ = 0.42 GeV cm^^; in- 
stead when po = Pq™ = 0.20 GeV cm""^, L may go up 
to 15 kpc with a range of S which progressively shrinks 
to (5 ~ 0.70 — 0.75, when L increases. 

We have also shown that future measurements of galac- 
tic antiprotons and antideuterons will offer, together with 
the upcoming data from DAMA/LIBRA, very interesting 
perspectives for further investigating the light neutralino 
populations selected by the annual modulation data. In 
case of models with a corotating halo or with triaxial spa- 
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FIG. 5: Antideuteron flux at D kinetic energy = 0.23 GeV/n, as a function of the WIMP mass and for a cored isothermal 
halo. Notations are as in Fig. O except for the horizontal lines, which here refer to estimated sensitivities to antideuterons of 
the GAPS (dashed) and AMS (dotted) detectors. 



tial distributions, not investigated in the present paper, 
also heavier neutralinos can be involved. 
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VI. APPENDIX: THE SUPERSYMMETRIC 
MODEL 

The supcrsymmctric scheme we employ in the present 
paper is the one described in Ref. an ef- 

fective MSSM scheme (cffMSSM) at the electroweak 
scale, with the following independent parameters: 
M2, ^, tan/3, m^, TOg, TO[, A and R = M1/M2. Notations 
are as follows: ^ is the Higgs mixing mass parameter, 
tan/3 the ratio of the two Higgs v.e.v.'s, m^i the mass 
of the CP-odd neutral Higgs boson, is a squark soft- 
mass common to all squarks, mj is a slepton soft-mass 
common to all sleptons, A is a common dimcnsionless 
trilinear parameter for the third family, = Ai = Arrig 
and Af = Amj (the trilinear parameters for the other 
families being set equal to zero). 

Since we are here interested in light neutralinos, we 
consider values of R lower than its GUT value: Rgut — 
0.5. For definiteness, we take R in the range: 0.005 - 0.5. 

In the present paper the numerical analyses are per- 
formed by a scanning of the supersymmetric param- 
eter space, with the following ranges of the MSSM 
parameters: 30 < tan/3 < 50, 100 GeV < < 



300 GeV, 100 GeV < M2 < 1000 GeV, 100 GeV < 
mq,mi < 1000 GeV, 90 GeV < tua < 150 GeV, -3 < 
A<3. 

The following experimental constraints are imposed: 
accelerators data on supersymmetric and Higgs boson 
searches (CERN e+e" coUider LEP2_|23 and Colhder 
Detectors DO and CDF at Fermilab_|28[ ) ; measurements 
of the b s + ^ decay process [29|: 2.89 < B{b — > 
5+7)- 10^^ < 4.21 is employed here: this interval is larger 
by 25% with respect to the experimental determination 
I29II in order to take into account theoretical uncertainties 



in the SUSY contributions [3l| to the branching ratio of 
the process (for the Standard Model calculation, we em- 
ploy the recent NNLO results from Ref. [s^]); the u ppe r 
bound on the branching ratio BR{B^ — > /i~ -|- [32]: 
wclakc BR{B1 ji^+ii^) < 1.2-10^^; measurements of 
the muon anomalous magnetic moment = (g^ — 2)/2: 
for the deviation Aa^ of the experimental world aver- 
age from the theoretical evaluation within the Standard 
Model we use here the range —98 < Aa^ • 10^^ < 565, 



derived from the latest experimental [32 
iil data. 
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